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ABSTRACT 

In order to characterise the cloud structures responsible for the formation of high-mass stars, we present Herschel observations of the DR21 
environment. Maps of the column density and dust temperature unveil the structure of the DR21 ridge and several connected filaments. The ridge 
has column densities larger than 10 23 cm -2 over a region of 2.3 pc 2 . It shows substructured column density profiles and branching into two major 
filaments in the north. The masses in the studied filaments range between 130 and 1400 M© whereas the mass in the ridge is 15000 M©. The 
accretion of these filaments onto the DR21 ridge, suggested by a previous molecular line study, could provide a continuous mass inflow to the 
ridge. In contrast to the striations seen in e.g., the Taurus region, these filaments are gravitationally unstable and form cores and protostars. These 
cores formed in the filaments potentially fall into the ridge. Both inflow and collisions of cores could be important to drive the observed high-mass 
star formation. The evolutionary gradient of star formation running from DR21 in the south to the northern branching is traced by decreasing dust 
temperature. This evolution and the ridge structure can be explained by two main filamentary components of the ridge that merged first in the 
south. 
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1. Introduction 

The properties of dense structures in interstellar molecular 
clouds, e.g., the density and temperature distribution, chemical 
composition, and dynamics, are expected to determine the de- 
tailed outcome of the star formation process that varies between 
an isolated low-mass star and an OB star cluster. While previous 
observation s have shown that clouds are filamentary (e.g., Bally 
et al.|1987|), recent Herschel studies have revealed the ubiq uity 



(e.g., |Miviile-Deschenes et al.|[20T0l |Henning et al.||2010| ) and 
key importance of filaments for star formation. Notably, grav- 
itationally un stable filaments are found to fragment into cores 
( [Andre et al.||2010| ), and filaments generally show narrow cen- 
tral widths of about 0.1 pc ( Arzoumanian et al.|2011| ). Resulting 
self-gravitating cores with sizes of 0.01-0.1 pc and masses up to 
10 M© are best candidates to form low- to possibly intermediate- 
mass stars (see rev iews by Di Fran cesco et al.|[2007 
Thompson etal.|2007 ) 



Ward- 



It is an important question which cloud structures lead to the 
formation of high-mass stars (>10M o ) almost exclusively found 
in stellar clusters or associations (see review by Zinnecker & 



|Yorke|2007|). The Herschel imaging survey of OB Young Stellar 



objects (HOBYS, Motte et aT1|2010| ) observes massive molec- 
ular cloud complexes within 3 kpc distance to probe the cloud 
environment of OB star- forming cores, their statistical evolu- 



tion, and the effects of feedback o n parental clouds. Studying 
the Vela C region, [Hill et al. ( 2011| ) suggested that "ridges", i.e., 
massive, gravitationally unstable filamentary structures of high 
column density (Nh 2 > 10 23 cm -2 ) that dominate their environ- 
ment could be preferential s ites of high-mass star formation (cf. 
Ngu yen Luong et aL|2011| ). The ridge in Vela C shows a com- 
plex substructure that could result from stellar feedback ( [Minier 
|et al.|2012] ). Further, intersecting filaments appear to mark sites 
of stellar cluster formation in the Rosette cloud ( Schn eider et al. 
[2012] ). 

This study focusses on the DR21 ridge (also called DR21 
filament), the densest and most massive cloud structure in the 
Cygnus X r egion ( [Schneider et al.|2006[|Motte et al.|2007[|Roy 
|et al.||2011| ) at a distance of 1.4 kpc ( [Rygl et al.||2012|). It host s 
the embedded HII region DR21 (e.g., |Roelfsema et al.|fT989]) 
the m aser source DR21(OH), an d mas sive protostars (Bontemps 
|et al.|[20T0| ). |Schneider et al.| ( [201 0| ) analysed molecular line 
emission and identified three "sub-filaments" F1-F3 that con- 
nect to the ridge. Velocity gradients suggest that material is 
transported along them towards the ridge, and in the case of F3 a 
bend and possible direct connection to the DR21(OH) clump is 
traced (cf. |Csengeri et al.|201 1| ). This behaviour supports the sce- 
nario that sustained accretion of inflowing material plays a role 
in building up massive clumps and cores and setting the stage 
for high-mass star formation, as suggested by numerical sim- 
ulations (e.g ., |Balsara et al.||2001 Banerjee et al.||2006 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA (Pilbratt et al. 2010). 



Smith 



et al. 2011| ). Here we exploit the unprecedented sensitivity of 



Herschel far-infrared and submillimeter continuum imaging to 
trace the detailed column density structure in this region. 
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Fig. 1: Herschel maps of the DR21 environment showing a) 70/mi emission, b) column density, and c) dust temperature. The DR21 
ridge is delimited roughly by the Nh 2 = 10 23 cm -2 contour plotted in panels b) and c). The filaments selected using DisPerSE (see 
Section [31 are named and marked with dots along their crests in b) and c). 



2. Herschel observations and data reduction 

Cygnus X North was observed in the parall el scan map mode 
with PACS (Po glitschetal.|2010) and SPIRE ( [Griffin etal.|20l0) 
on Dec 18, 2010. To diminish scanning artifacts, two nearly 
perpendicular coverages of 2.8x2.8° were obtained in five pho- 
tometric bands at 70, 160, 250, 350, and 500 yum with a scan 
speed of 20 // /sec. The data were reduced using the Herschel 
Interactive Processing Environment and the Scanamorphos soft- 
ware (Roussel 2012, for creating PACS maps). 

The 70 fim map traces in particular heated dust towards star- 
forming sites (Fig. [T]a). To quantitatively estimate the distribu- 
tion of the cold dust which generally represents most of the dust 
mass, we created maps of column density Nh 2 and dust tem- 
perature T d in the way described in Hill et al. ( 201 1| [2QT2] ), in- 
cluding offsets to rec over the absolute intensity level following 
Bernard et aL] ( [20T0| ). Computed from the 160, 250, and 350 jim 
bands, the maps provide an angular resolution of 25"(0.17pc at 
1.4 kpc). They have undefined values for pixel groups towards 
DR21, DR21(OH) (white pixels in Fig. [T]b), and W75N (see 
Fig.|B.l|) due to saturation at 250 and 350 yum. 



3. The DR21 filaments 

Shown in Fig.[T]a, compact 70 jim sources that are protostar can- 
didates cluster along the DR21 ridge oriented north-south, with 
the most prominent peak being DR21 itself. Supposedly lying 
at the same distance, they show a strong northward decrease in 
luminosity. Many filamentary streamers from the ridge to e.g., 
the north-west and west, are present, most of which were also 
detected in the m id-infrared with Spitzer (Mar ston et al.||2004} 
Hor a et al.|20 09). They mainly correspond to low column den- 
sity structures (Fig. [T]b), but several filaments are prominent. 
The extent of the DR21 ridge can be roughly defined by the 
NH 2 = 10 23 cm -2 conto ur enclosing an area o f 2.3 pc 2 . Notably, 
the "sub-filaments" of |Schneider"etaLl ( [20T0] ) labeled Fl and F3 
are recovered, both being resolved into nearly parallel northern 
and southern components that join close to the ridge. The north- 
ern part of the DR21 ridge shows two extensions in column den- 
sity ("rabbit ears") to the north and north-west (Fl) which coin- 
cide with the coldest regions where the dust temperature drops 



as low as 14 K (Fig. [T]c). The filaments are visible in the dust 
temperature map as structures of lower central temperature rel- 
ative to the background of -19 K. The background column den- 



sity level is ~10 cm with a standard deviation of ~10 cm 
mainly due to cirrus structure. 

In the first step to characterise the DR21 filaments, we ap- 
plied the DisPerSE software (Sousbie||2011|) on the column 



density map to identify the filament crests (cf. |Arzoumanian 



|et al.||2011| Hill et aL||2011| ) using a persistence threshold of 
3 x 10 2i cm . Unlike for these works, possible line-of-sight con- 
fusion needed to be addressed in our analysis due to the low- 
density foreground cloud Cygnus Rift and the close-by W75N 
cloud component ( e.g., Schneider et al.||2006|). Also the str ong 
DR21 outflow (e.g., |Richardson et al.|1986[|White et al.|2010) is 



avoided. Subsequent studies will analyse the combined column 
density, temperature, and kinematical data for a larger region, 
but here we conservatively limited this study to filaments that 
(1) have main velocities in 13 CO coherent with the ridge and no 
strong secondary velocity component^] and (2) have crest col- 
umn densities above 1.5 x 10 22 cm~ 2 ; the resulting selection is 
shown in Fig.[T] At positions along their crests spaced by half the 
beam FWHM we extracted perpendicular column density pro- 
files extending over 2.5 pc in each direction, wide enough so that 
the outer profiles are sufficiently flat and trace the background. 
For background column densities, we adopted the profile min- 
ima that were median- smoothed over 2.5 beam FWHM to lower 
the influence of neighbouring structures. 

Crest column density and dust temperature are anti- 
correlated as anticipated for all filaments with the possible ex- 
ception of SW (Fig. [2]a, b). Contrary to the others, the lat- 
ter shows slightly decreasing central temperature outward from 
the DR21 ridge. This filament also shows the weakest outward 
column density decline. In general, lower temperatures are ex- 
pected for structures of higher column density due to the in- 
creased shielding from the interstellar radiation field. However, 
the SW filament is running towards the warmer area around 
DR21 (Fig. [T]a, c) and possibly into the outflow cavity. Its in- 



a We used the |Schneider et al.| ( |2010| > FCRAO 13 CO(1-0) data and 
excluded filament segments showing additional velocity components 
above the 20% level of the primary peak. 
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Table I: Average properties of the DR21 filaments and ridge 





S 


sw 


F3S 


F3N 


F1S 


FIN 


N 


Ridge 


Length (pc) 


2.6 


1.6 


1.3 


2.7 


1.3 


6.9 


1.1 


4.1 


Mean crest column density a (10 22 cm -2 ) 


4.4 


1.2 


1.6 


2.1 


7.7 


1.9 


4.8 


41.6 


Integration radius (pc) 


0.58 


0.41 


0.24 


0.24 


0.48 


0.39 


0.49 


1.15 


Total mass (M©) 


1350 


210 


130 


400 


960 


1210 


610 


15210 


Mean line mass (M©/pc) 


500 


120 


100 


140 


700 


170 


520 


3670 


Weighted mean central width^ (pc) 


0.28±0.07 


0.31+0.07 


0.27+0.08 


0.32+0.11 


0.26+0.07 


0.33+0.16 


0.34+0.17 


0.34+0.12 



Notes. (a) Background- subtracted column density. (b) The fit error in width <x w was used for weighting data points with 1/cr 2 



verted density-temperature relation could thus be an effect of 
heating by DR21. 

The perpendicular filament column density profiles gener- 
ally show a flattened central part and non-Gaussian, roughly 
power-law wings. To examine the central width we performed 
a weighted fit of a Gaussian curve plus a constant level (cf. 
|Arzo umanian et al. 2011 ). As weights a Gaussian curve with 
0.17 pc FWHM was used to ensure the reproduction of the cen- 
tral flattening. The central width of each profile was derived as 
deconvolved fitted Gaussian FWHM. The weighted mean val- 
ues of the profile central widths lie between 0.26 and 0.34 pc 
(Table [T}. These values are beyond the rang e around the typical 
~0.1 pc found by Arzoumani an et al. (2011), probably only par- 
tially caused by the lower spatial resolution. However, no corre- 
lation between central width and crest column density is present, 
neither for profiles nor for the mean values. Better statistics in 
the high-column density regime is needed to clarify a possible 
dependency. 

To constrain the masses of the filaments we estimated an 
outer radius for each profile integration. A threshold per pro- 
file was derived by building the cumulative integral over radial 
distance from the crest and determining the first minimum of its 
derivative. For each filament we adopted its median threshold as 
integration radius, ranging from 0.24 to 0.58 pc (Table [T}. 

In theory, the thermal stability of filaments is to fi rst order de- 
termined by the critical mass per unit lengtrj^(e.g., Inutsuka & 
Miyama 1992]). Observations indica te that this is a proxy for the 
general stability (Andre et aL 2010| ). Most filament segments are 
thermally supercritical (Fig. 2]c). The outer segments (beyond 
1 pc) of the SW, F3N, and F3S filaments are only supercritical 
to within a factor of 3, and the longest, FIN filament shows a 
marginally subcritical segment. For thermally supercritical fil- 
aments, fragmentation and local spherical collapse occur faster 
than g lobal collapse in the absence of o ther stabilising mecha- 
nisms ( |Pon et al.|201 l[|Tbala et al.|2012| ). In this case we thus ex- 
pect the DR21 filaments to form cores and protostars. The S fila- 
ment shows several local maxima of crest column density seper- 
ated by about 1 pc (Fig.|2]a). Less regular and prominent, several 
peaks are also seen for the other filaments. At 70/mi there is 
one extended source towards the S filament (Figs. [T]a, |A.l| a). 
Compact emission from protostar candidates is present towards 
the junction point of F3S and F3N, and towards the FIN and 
N filaments. The 250 yum map shows additional compact star- 
less/prestellar core candidat es towards a ll filaments, one of them 
identified as dense core by |Motte et al. ( 2007| ) (see Apps.[A|[B]). 
This confirms that core and star formation is ongoing within the 
thermally supercritical filaments. 




b M^ e = 2c 2 s /G with the isothermal sound speed c s = ^jKYjJi^ 
and T d at the crest position, resulting in Mjf" ^30M o /pc. 



1.0 1.5 2.0 2.5 3.0 4.0 5.0 6.0 7.0 
Distance along subfilament (pc) 

Fig. 2: DR21 filament properties over distance along their crest 
from the DR21 ridge. The distance scale is compressed be- 
yond 3 pc. a) Crest column density (not background- subtracted), 
dashed curves show the background column density, b) Crest 
dust temperature, c) Mass per unit length, dashed lines gives the 
critical value. 



4. The DR21 ridge - an intersection of filaments 

The DR21 ridge hosts 22 dense cores and -1/4 of the mas- 
sive ones in the |Motte et aL] ( |2007| ) sample of the whole cloud 
complex, and thus represents the probable region to host the 
most massive forming star cluster in Cygnus X. We used their 
MAMBO 1.2 mm map to interpolate the missing pixels in the 
column density map and derived the mean ridge properties us- 
ing a crest of 4pc length following its main peaks (Fig. [3]a, 
Table [TJ. If it were regarded as an individual filament, it would 
be highly un stable, in accord with th e global collapse signature 
observed by Schneider et al. (2010) and the dense core detec- 
tions (Fig.[3]a). 

The chain of cores dominates the structure of the ridge ex- 
cept between DR21 and DR21(OH), where a crest column den- 
sity dip is present (location c in Fig. 3|. Less clear than the north- 
ern branching, a possible additional filamentary extension coin- 
ciding with a core branches off to the south-east at the south- 
ern end of the ridge (at the very bottom of Fig. [3]a). As illus- 
trated in Fig.[3]b for location b, corresponding to the dense core 
N40, the column density profile of the ridge generally shows 
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^ Averaged profile 
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L 
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Radius (pc) 



mation involving a core merg er process was also suggested by 
observations of NGC 2264-C jPeretto et al.|2Q06| ). 

The Herschel observations emphasize the evolutionary gra- 
dient along the ridge: Beyond DR21, the 70 fim luminosity of 
protostars strongly decreases northward, and the dust tempera- 
ture shows lowest temperatures towards the northern part. The 
substructure of the DR21 ridge suggests that it was formed by 
the merging of individual narrow, intersecting filaments. The 
pronounced elongation and branching into filaments Fl and N 
with the highest mass per unit length indicate that the major 
components may be roughly parallel filaments oriented north- 
south. At present, the merging could have advanced furthest in 
the southern part and less in the north, where the two compo- 
nents appear separated. Extrapolating this scenario, we expect 
that the northern filaments will lead to the assembly of one or 
more further massive clumps. This study suggests that high-mass 
star-forming ridges could be second-generation cloud structures 
formed via dynamical merging of gravitationally unstable fila- 
ments. 



Fig. 3: a) Column density map (in cm -2 ) of the DR21 ridge with 
the crest as black dots and the black Nh 2 = 10 23 cm -2 contour. References 
Note that white pixels are undefined due to saturation. Blue el- 
lipses give the po sition and extent (FWHM = semiaxes) of the 
Motte et aL| ( |2007| ) dense cores, b), c) Radial column density pro- 
files towards the crest locations b and c marked in yellow/blue in 
a). The radial profiles are averaged over the marked crest dots. 
The yellow areas show the ±1 cr ranges. 



narrow peaks and smoothly decreasing wings towards the cores. 
However, the profile towards location c shows at least one ad- 
ditional peak of 9xl0 22 cm -2 at 0.4 pc (Fig.[3]c), indicating that 
there the ridge consists of more than one individual filament. 
Considered together with the branching of the northern ridge 
into the "rabbit ears" (the Fl and N filaments) and the possible 
southern branching, this finding suggests that the DR21 ridge is a 
complex intersection of several individual filaments. The domi- 
nating filament along the crest shows a central width of -0.34 pc, 
not significantly larger than the central widths of the filaments 
(Table[T]). 



5. Star formation in the DR21 ridge 

Schneider et al.|([20T0|) present velocity gradients that suggest ac- 



cretion of the filaments onto the DR21 ridge. Dust temperature 
dips are seen towards the filament-ridge connections (Fig. [2]b), 
indicating that inflowing material has cooled at these locations. 
At present, the total mass in the selected filaments is about 1/3 of 
the ridge mass. It thus appears that the mass assembly process of 
the DR21 ridge is in a late stage, possibly driven by the gravita- 
tional potential of the ridge itself, and specifically the build-up of 
massive cores in the DR21/DR21(OH) clumps could have been 
supported by continuous accretion from the S/SW and F3S/F3N 
filaments, respectively. The filaments are for the most part grav- 
itationally unstable and forming cores and prot ostars, in contrast 
to the striations and "sub-filaments" in Taurus (Golds mith et al.| 
2008, Palmeirim et al. in prep.) or Aquila ( Andr e et al.||2010| ). 
The high masses of the cores in the ridge could thus not only be 
due to the filament flows, but also due to the merging with frag- 
ment cores of the filaments that form dense, small-scale conver- 
gent flows. For DR21(OH), Csen geri et al.| ( |20lT] ) showed possi- 
ble fragmentation of the inflowing material. High-mass star for- 
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Appendix A: The S and F1N filaments in detail 

Figure [ATT] shows a zoom on the S and FIN filaments. Several 
compact, more or less elongated cores are seen in the 250 /mi 
and column density map. The marked dense core N23 has an 
estimated mass of 13 M . 
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Fig. A.l: Detailed view of filament S (panels a, b, c) and FIN (d, e, f). The 70 and 250 /mi units are MJy/sr, column density unit 
is cm" 2 . The black cross in panel a) indicates an extended 70 fim source. The blue ellipse in panels d), e), f) gives the position and 



extent (FWHM = 



: semiaxes) of the dense core N23 of |Motte et al. (2007 ). 



6 



M. Hennemann, F. Motte, N. Schneider et al.: The spine of the swan, Online Material p 1 



Appendix B: Herschel maps of DR21 and Cygnus X 
North 

Figure |BJ"] shows the DR21 environment in the Herschel maps 
of the HOBYS Cygnus X North observations together with a 
RGB colour composite image. 
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Fig. B.l: Herschel maps of the DR21 environment in Cygnus X North obtained by HOBYS. Map units are MJy/sr. The last panel 
shows a RGB colour composite image using the SPIRE 250 yum (red), PACS 160yum (green), and PACS 70/mi (blue) maps. 



